Search for d* Dibaryon by Double-radiative 
Capture on Pionic Deuterium 



P. A. Zoinierczuk a , T. P. Gorringe a , M. D. Hasinoff b , 
M. A. Kovash a , S. Tripathi a , D. H. Wright ' 1 

^University of Kentucky, Lexington, KY 40506, USA 
b University of British Columbia, Vancouver, B.C., Canada V6T 1Z1 
c TRIUMF, 4004 Wesbrook Mall, Vancouver, B.C., Canada V6T 2A3 

Abstract 

We report a search for d* dibaryon production by double-radiative capture on pionic 
deuterium. The experiment was conducted at the TRIUMF cyclotron using the 
RMC cylindrical pair spectrometer, and detected 7-ray coincidences following pion 
stops in liquid deuterium. We found no evidence for narrow dibaryons, and obtained 
a branching ratio upper limit, BR < 6.7 x 10~ 6 (90% C.L.), for narrow d* production 
in the mass range from 1920 to 1980 MeV. 
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1 Introduction 



At present the deuteron is the only established particle with a baryon num- 
ber B = 2. However a large number of theoretical predictions for additional 
dibaryons have been published in the literature. The predictions have exploited 
both quark-gluon and hadron viewpoints, and include objects such as exotic 
multi-quark states and molecular-like baryonic states. Obviously the experi- 
mental discovery of another dibaryon would provide new insight into hadron 
and quark-gluon dynamics at the GeV scale. 

A recent claim for dibaryon production has been published by the Di27 collab- 
oration at the JINR phasotron [1]. The collaboration measured the proton- 
proton double bremsstrahlung reaction pp — > ppjj by directing 216 MeV 
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protons onto a liquid hydrogen target and recording 7-ray coincidences in 
two Csl/Nal detector arrays. The authors observed an intriguing structure in 
the background-subtracted energy spectrum of the 7-ray coincidence data. It 
comprised a relatively narrow peak centered at about 24 MeV and a relatively 
broad peak centered at about 60 MeV. They attributed the structure to the 
d* dibaryon with mass 1956 ±6 MeV and width ^ 8 MeV. They hypothesized 
that d* dibaryons were first produced via the two-body process pp — > d*-y and 
then decayed via the three-body process d* — > ppj, thus contributing to pp 
double bremsstrahlung. 

The d* dibaryon of Khrykin et al. [1] would have an electric charge Q = +2, 
an isospin component T z = +1, and therefore an isospin of T ^ 1. The narrow 
width of the claimed d* state implies that its quantum numbers either forbid 
or strongly suppress the NN decay mode. Thus, if the d* is isovector (T = 1), 
it must have spin-parities J 71 " = l + ,3 + ,5 + etc. Moreover Khrykin et al. have 
argued that (J n ,T) = (1 + , 1) is the most natural choice for the d*'s quantum 
numbers, being the lowest spin-isospin pp-decoupled state with zero orbital 
angular momentum. 

Recently Gerasimov [2,3] suggested that double-radiative capture on pionic 
deuterium 

n~d — > 717177 (1) 



is an excellent candidate for further investigations of the dibaryon's existence. 
In Eqn. 1 the d* (T z = — 1) dibaryon is first produced via radiative capture 
iT~d — > <i*7 and then disintegrates via radiative decay d* — > nn r ). Using a sim- 
ple model Gerasimov estimatedP] that the branching ratio for the <i*-mediated 
process might be as large as 0.5%. This yield would exceed by 100 times the ex- 
pected two-photon branching ratio for non-resonant double-radiative capture 
in pionic deuterium (see [2,4]). 

A number of 7-ray experiments on pionic deuterium have already been con- 
ducted, so is it possible for the signatures of the d* dibaryon to have been 
missed? Singles 7-ray data on pionic deuterium are available from the TRI- 
UMF experiments of Highland et al. [5] and Stanislaus et al. [6] and the PSI 
experiment of Gabioud et al. [7]. Unfortunately, because of the large branch- 
ing ratio BR = 0.26 for the single radiative capture reaction, which yields a 
~130 MeV 7-ray peak with a large low-energy tail, a small contribution from 
dibaryon production with a branching ratio ^10~ 2 is not excluded (the exact 
sensitivity of the singles 7-ray experiments is a strong function of the mass 
and the width of the dibaryon) . Coincidence 7-ray data on pionic deuterium 

2 Gerasimov assumed a NA bound state for the d* dibaryon and a Kroll-Ruderman 
graph ttN — > 7A for the production mechanism. 
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are also available from the TRIUMF experiment of MacDonald et al. [8] . How- 
ever this experiment was designed for back-to-back photons originating from 
7r° — > 77 decay following 2 H(7r~, ir°) charge exchange. Consequently their sen- 
sitivity to two-photon events from d* dibaryon production in double-radiative 
capture was very low. 



2 Experimental Setup 

Herein we report a dedicated search for dibaryon production via double- 
radiative capture on pionic deuterium. The experiment was conducted on the 
M9A beamline at the TRIUMF cyclotron using the RMC pair spectrometer 
(see Fig. 1). 

The beamline delivered a negative pion flux of 5 x 10 5 s" 1 with a central 
momentum of 81.5 MeV/c and a momentum bite Ap/p of about 10%. The 
beam was unseparatedQ having a e/n ratio of about 14/1 and a /a/it ratio of 
about 1/1. The incoming pions were counted in a 4-element plastic scintillator 
telescope and stopped in a 2.5 liter liquid deuterium target. 

The outgoing photons were detected by electron-positron pair production in 
a 1 mm cylindrical lead converter and e + , e~ tracking in a cylindrical drift 
chambers. A 1.2 kG axial magnetic field was used for momentum analysis 
and concentric plastic scintillator rings were used for fast triggering. The trig- 
ger scintillator package consisted of the A-ring (just inside the Pb converter 
radius), the C-ring (just inside the multiwire drift chamber radius), and the 
D-ring (just outside the drift chamber radius). For more information on the 
spectrometer see Wright et al. [9], note that in this experiment we moved 
the Pb converter from just inside the C-counter radius to just outside the 
A-counter radius. 

The two-photon trigger was based upon the multiplicities and topologies of 
hits in the trigger scintillators and the drift chamber cells. The scintillator 
trigger required zero hits in the A-ring, two or more hits in the C-ring, three 
or more hits in the D-ring, and a C-D topology consistent with the conversion 
of two 7-rays. The drift chamber trigger imposed minimum values for number 
of the drift cell hits or the drift cell clusters in each drift chamber layer. 

During a four week running period we accumulated 7-ray coincidence spectra 
from a total of about 3.8 x 10 11 pion stops in the deuterium target. We also 
collected data from it~ stops in liquid H 2 for setup and calibration. 
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3 Data Reduction 



The significant backgrounds involved (i) true coincidences from tt° 77 de- 
cay following (7r~,7r°) charge exchange on X H contamination, (ii) accidental 
coincidences from two ir~d — > rmq events following two pion stops in one 
beam pulse (hereafter denoted it-it events) Q and (iii) accidental coincidences 
between a delayed radiative \i decay 7-ray and a prompt radiative it cap- 
ture 7-ray (hereafter denoted fi-ir events). The tt° decay background yielded 
photon-pairs with opening angles of cos 6*12 ^ —0.76 and summed energies 
of E sum = E y i + E^ 2 — rn n . The accidental coincidence backgrounds yielded 
photon-pairs with opening angles of —1.0 ^ cos 9 12 ^ +1.0 and summed 
energies of up to ~ 180 MeV for fi-ir accidentals and ~ 260 MeV for it-it 
accidentals. 

In analyzing the data a number of cuts were applied to select the two-photon 
events. A tracking cut imposed minimum values for the number of points in 
the tracks in the drift chamber and maximum values for the chi-squared of fits 
to the tracks. A photon cut required that the electron-positron pairs intersect 
at the lead converter and that the photon pairs originate from the deuterium 
target. The total number of photon pairs, i.e. events surviving the tracking 
cuts and photon cuts, was 2.3 x 10 5 . These photon pairs (shown in Fig. 2) 
are dominated by the real 7-7 coincidences from tt° decays and accidental 7-7 
coincidences from it-it events and \x-tt events. The accidental 7-7 background 
is clearly seen in the summed energy spectrum as events with E > 150 MeV. 
The 7T° decay background is clearly seen in the opening angle spectrum as 
events with cos #12 < —0.6. In order to remove these backgrounds, a beam 
counter amplitude cut was applied to reject the tt-tt accidentals and a en- 
counter timing cut to reject the accidentals. Finally an opening angle cut 
was used to reject the background from ir° decay. 

A total of 370 two-photon events were found to survive the beam counter 
amplitude cut, C-counter timing cut, and photon opening angle cut. Q Their 
photon opening angle and individual energy spectra are shown in Fig. 3. The 
opening-angle spectrum is a broad continuum and shows the opening angle 
cut of cos 9 12 > —0.2 we employed to remove the events from tt° decay. The 
individual energy spectrum is a broad continuum with a low-energy cut-off at 

4 The pion beam had a micro-structure with a pulse width of 2-4 ns and a pulse 
separation of 43 ns. For an incident flux of 5 x 10 5 s _1 the probability for more than 
one pion arriving in a single beam pulse was about 1.1%. 

5 Actually a small contribution of 22 ± 5 tt° decay background events and 24 ± 2 
accidental coincidence background events is present in Fig. 3. The residual tt° decay 
background was estimated using the number of tt° events with cos #12 < —0.45 and 
the residual accidental coincidence background was estimated using the number of 
tt-tt / fi-ir events with E sum > 150 MeV. 
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about 25 MeV due to the acceptance of the spectrometer. 

The 370 two-photon events in Fig. 3 are entirely consistent with non-resonant 
capture in pionic deuterium. Specifically, using the measured branching ratio 
for non-resonant double radiative capture on hydrogen of (3.05±0.27±0.31) x 
10~ 5 [10], and naively assuming the ratio of single radiative capture to double 
radiative capture to be identical on a hydrogen target and a deuterium target, 
we would expect about 580 ± 110 events from non-resonant n~d — > 777m. 
In addition the measured angle and energy distributions are in good agree- 
ment with our Monte Carlo simulations of the non-resonant process using the 
theoretical model of Beder [11]. The expected signature of dibaryon events, a 
monoenergetic peak from the production process ir~d — > d*^ and a three-body 
continuum from the decay process d* — > 71717, is not seen. 



4 Dibaryon Sensitivity 

To determine the detection efficiency for dibaryon events we used a Monte 
Carlo computer program. The program incorporated the detailed geometry 
of the RMC detector and detailed interactions of the various particles. Our 
program was based on the CERN GEANT3 package [12] and is described in 
more detail in Ref. [13]. 

The simulation was tested by measurements of the detector response with a 
liquid H 2 target. Negative pion stops in hydrogen provide a well known source 
of photon pairs from (i) ir~p — > 7T°n charge exchange followed by 7r° — > 77 
decay and (ii) accidental 7-7 coincidences from multiple tt stops. We found 
the energy-angle distributions from experiment and simulation to be in good 
agreement. The absolute detection efficiencies from experiment and simulation 
differed by 4% (the run-to-run variation of the spectrometer acceptance was 
< ±6%). This difference was attributed to detector inefficiencies that were 
present in the measurement but were absent in the simulation. A multiplicative 
correction factor F = 0.96 ± 0.06 was therefore incorporated to account for 
this. 

As described earlier, the d* signature is a monoenergetic 7-ray peak from 
dibaryon production, ir~d — > c?*7, and a three-body continuum from dibaryon 
decay, d* — > nny. The production 7-ray energy is approximately E ~ M^-a — 
M d * and increases from E ~ 35 MeV for M d * = 1980 MeV to E ~ 90 MeV 
for M d * = 1920 MeV. The decay 7-ray spectrum is peaked near the end-point 
energy, E ~ Md* — 2m n , as typically the two neutrons carry only a small 
fraction of the available energy. However the detailed shape of the three-body 
spectrum is dependent on the nn-final state interaction and consequently on 
the d* quantum numbers. 
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Representative simulations of dibaryon signatures for Md* = 1920 and 1956 
MeV, assuming a 3-body phase space distribution for d* — * nn / ~f decay, are 
shown in Fig. 4. Note that for M d * = 1920 MeV the production and decay 
7-ray spectra are separated whereas for M d * = 1956 MeV the production and 
decay 7-ray spectra overlap. Because the exact line shape of the decay 7-ray 
is not known, we obtained our limits on d* production by determining limits 
on the production 7-ray yield. Specifically we fit the sum of a polynomial 
function, parameterizing the non-resonant background, and a Gaussian peak, 
accounting for d* production, to the 7-7 coincidence energy spectrum in Fig. 3. 
The Gaussian peak centroid was stepped from E = 35 to 90 MeV and the 
Gaussian peak width was fixed at the instrumental resolution of the RMC 
detector (i.e. a =4-5 MeV for E =40-80 MeV). This procedure yielded a limit 
on the number of dibaryon events as a function of dibaryon mass. We then 
converted the event limit on dibaryon production to a branching ratio limit 
on dibaryon production via 

B - K = N^-eAQ-F-C (2) 

where Nj* is the limit on the dibaryon events, N n - is the number of live-time 
corrected pion stops, and eAQ ■ F is the detector acceptance. The appropri- 
ate acceptance was obtained using the Monte Carlo for dibaryon production 
and assuming a three-body phase space distribution for d* — > nn / -f decay 
(cf. Fig. 5). As discussed earlier the factor F accounts for detector inefficien- 
cies which are present in the experiment but are absent in the simulation. The 
factor C = 0.85 ±0.01 accounts for the fraction of incident pions that stopped 
in deuterium (see Wright et al. [13] for details). 

The resulting branching ratio upper limit on d* production in 7r~d capture 
was found to be BR < 6.7 x 10~ 6 (90% C.L.) for rf*'s in the mass range 
of 1920 to 1980 MeV and width of < 10 MeV (see Fig. 5). In particular, 
we observed no evidence for a narrow dibaryon of mass M = 1956 MeV as 
claimed by Khrykin et al. [1]. However above and below this mass range, our 
experimental sensitivity rapidly deteriorates due to the energy cut-off in the 
spectrometer acceptance. 



5 Summary 

In summary we have found no evidence for narrow dibaryon production in 7i~d 
capture. Our upper limit on dibaryon production, BR < 6.7 x 10~ 6 (90% C.L.), 
is several orders of magnitude below the yield estimate of Gerasimov [2]. Our 
null result is consistent with the null result of Calen et al. [14] obtained from 
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pp bremsstrahlung measurements using the WASA detector at the CELSIUS 
storage ring. 

We wish to thank the staff of the TRIUMF laboratory for their support of 
this work. In particular we acknowledge the help of Dr. Renee Poutissou on 
the data acquisition and Dr. Dennis Healey on the cryogenic targets. In ad- 
dition we thank Prof. Sergo B. Gerasimov for helpful and stimulating discus- 
sions, Prof. David S. Armstrong for carefully reading the manuscript and the 
National Science Foundation (United States) and the Natural Sciences and 
Engineering Research Council (Canada) for financial support. 



References 



[1] A. S. Khrykin et al, Phys. Rev. C 64, 4002 (2001). 

[2] S. B. Gerasimov, Proc. of the Int'l. Workshop "Hadronic Atoms and 

Positronium in the Standard Model", Dubna, 26-31 May 1998, JINR-E2-98- 
254, Dubna 1998 p.91; arXiv.org jnucl-th /98080"7C . 

[3] S. B. Gerasimov, Proc. of the XV Int'l. Seminar on High Energy Physics 
Problems, Dubna, September 25-29, 2000, JINR-E1, 2-2001-291, Dubna, 2001, 
Vol. II, p.166. 

[4] P. A. Zolnierczuk, Acta Phys. Pol. B 31 (2000) 2349. 

[5] V. L. Highland et al, Nucl. Phys. A 365, 333 (1981). 

[6] S. Stanislaus et al, Phys. Lett. B 219, 237 (1989). 

[7] B. Gabioud et al, Nucl. Phys. A 420, 496 (1984). 

[8] R. MacDonald et al, Phys. Rev. Lett. 38, 746 (1977). 

[9] D.H. Wright et al, Nucl. Instrum. Methods A320, 249 (1992). 

[10] S. Tripathi et al, D. S. Armstrong, M. E. Christy, J. Clark, T. P. Gorringe, 
M. D. Hasinoff, M. Kovash, D. H. Wright and P. A. Zolnierczuk, nucl- 



| ex/0204013| , Phys. Rev. Lett., accepted for publication. 
[11] D. Beder, Nucl. Phys. B156, 482 (1979). 

[12] R. Brun et al, CERN report no. DD/EE/84-1 (unpublished). 
[13] D.H. Wright et al, Phys. Rev. C 57, 373 (1998). 
[14] H. Calen et al, Phys. Lett. B 427, 248 (1998). 



7 




Fig. 1. The RMC spectrometer showing the deuterium target, lead converter, cylin- 
drical drift chamber, trigger scintillators and spectrometer magnet. 
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Fig. 2. The photon opening angle cos 9±2 (top) and summed photon energy spectra 
(bottom) for events passing the tracking cuts and photon cuts. The events are 
mainly real 7-7 coincidences from ir° decay and accidental 7-7 coincidences from 
7r-7r and fi-n events. 
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Fig. 3. The photon opening angle cos 9 12 (top) and individual photon energy spectra 
(bottom) for events passing all cuts. Our experimental data are denoted by full cir- 
cles and the dashed lines represent Monte Carlo simulations assuming non-resonant 
pion double-radiative capture on deuterium. 




Fig. 4. Examples of the individual 7-ray energy spectra for dibaryons with 
M = 1956 MeV (solid line) and M = 1920 MeV (dashed line) normalized using 
Eqn. 2 and the corresponding branching ratio upper limit taken from Fig. 5. We 
assumed a three-body phase space distribution for the d* — > jnn decay 7-ray spec- 
trum. 
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Fig. 5. The 90% C.L. dibaryon branching ratio upper limit versus the d* mass (full 
squares and left-hand scale) and the Monte Carlo acceptance versus the d* mass 
(dashed line and right-hand scale). 
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